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Abstract 



The reduced tension ad of the interface between the confined and the deconfined phase 
of SU (3) pure gauge theory is related to the finite size effects of the first transfer matrix 
eigenvalues. A lattice simulation of the transfer matrix spectrum at the critical temperature 
Tc = 1/Lt yields acd = 0.139(4)T^^ for Lj = 2. We found numerical evidence that the 
deconfined-deconfined domain walls are completely wet by the confined phase, and that the 
confined-deconfined interfaces are rough. 
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1 Introduction 



At high temperatures QCD does not confine color and chiral symmetry is not spontaneously 
broken. In the early universe, at about 10~^ seconds after the big bang, quarks and gluons 
have passed a transition to the low temperature confined phase in which chiral symmetry 
is spontaneously broken. Depending on the quark masses the transition may be a first or 
second order phase transition or just a cross over. At a first order phase transition the 
nucleation of bubbles of the confined phase from the supercooled high temperature quark- 
gluon plasma produces spatial inhomogeneities in the baryon density. If such inhomogeneities 
are large enough, they may infiuence the primordial nucleosynthesis of light elements which is 
completed within the first minute after the big bang. Small fiuctuations, on the other hand, 
are washed out and the baryon density at the time of nucleosynthesis becomes homogeneous. 
The size of the inhomogeneities is governed by the value of the interface tension between 
the high and the low temperature phase. It is convenient to introduce the reduced interface 
tension 

F 

(^cd=J^, (1-1) 

where F is the interface free energy, A is the interface area and T is the temperature. 
Computing is a nonperturbative problem which is numerically very difficult when quarks 
are present. Light quarks have a tendency to weaken the phase transition. To obtain an 
upper limit for the interface tension it is therefore sufficient to neglect quarks and to restrict 
oneself to a pure SU{3) gauge theory of only gluons. Then, as numerical lattice simulations 
have shown the phase transition is first order. The value of the reduced interface tension 
has been determined numerically at the critical temperature Tc = 1/Lt where Lt is the lattice 
extent in the euclidean time direction. For Lj = 2 the Boston group obtained acd = 0.12{2)T^ 
[0] while the Helsinki group quotes acd = 0.08 (2)T^^ 0. In this paper, using a completely 
different method, we obtain acd = 0.139(4)T^^. Closer to the continuum limit for = 4 the 
Boston group quotes acd = 0.027(4)T^^ [|[ . Taking this value as an estimate of the interface 
tension in real QCD with quarks the typical scale of inhomogeneities (e.g. the distance 
between nucleation centers) is at most a few centimeters 0. The proton diffusion length, 
on the other hand, is much larger (about 50 cm) such that the fiuctuations are washed out 
before primordial nucleosynthesis takes place. 

In previous numerical studies of the interface tension coexistence of the confined and 
the deconfined phase was enforced by keeping different parts of the lattice at different tem- 
peratures or by applying an external field. This breaks translation invariance and pins the 
interface at a certain position. The properties of a pinned interface will in general be different 
from the ones of the free interface one is interested in. To extract the interface tension one 
should first perform the infinite volume limit and then turn off the temperature gradient or 
the external field. In practice this is difficult because in a numerical simulation the lattice 
size is necessarily limited. This causes finite size effects which must be understood before one 
can extrapolate results reliably to the infinite volume limit. In this work we use a different 
strategy. We make use of the fact that the finite volume alone supports coexisting bulk 
phases when a cylindrical geometry is used. The lattices we use have two rather short x- 
and y-directions and one much longer 2;-direction. Then, without artificial temperature gra- 
dients or external fields, several bulk phases, aligned along the z-direction, coexist with each 
other, separated by confined-deconfined interfaces spanned in the short x- and y-directions. 
The presence of the interfaces infiuences the spectrum of the transfer matrix in the long 
z-direction. In particular, the lowest energy levels show a characteristic dependence on the 
interface area and on the reduced interface tension acd- We will use this finite size effect to 



extract the value of acd- This method has been appUed before to the four- dimensional ^ 
and to the three-dimensional Ising model 0. It has the advantages that no pinning of the 
interfaces is necessary and that finite size effects are well understood such that the results 
of numerical simulations can be extrapolated reliably to the infinite volume limit. 

For the pure glue theory, however, additional phenomena arise. Since in the high tem- 
perature phase the Z(3) center symmetry of the nonabelian gauge group is spontaneously 
broken there are actually three deconfined phases. At temperatures above Tc these phases 
coexist. They are separated by deconfined-deconfined domain walls with another reduced 
interface tension a^diT) which depends on the temperature. At very high temperatures 
(T ^ ^qcd) the reduced interface tension has been computed semiclassically p 

^UT) = (1.2) 

where g is the gauge coupling renormalized at the scale T. Frei and Patkos ^ have suggested 
that 

(^ddiTc) = 2acd, (1-3) 

i.e. at the transition temperature a deconfined-deconfined domain wall costs the same free 
energy as two confined-deconfined interfaces. In this situation — called complete wetting — 
the phase transition proceeds via domain wall splitting. Already slightly above the transition 
temperature the deconfined-deconfined domain walls split into two confined-deconfined in- 
terfaces by creating a complete wetting layer of confined phase between two bulk deconfined 
phases |]12|. Complete wetting is a critical phenomenon of interfaces which arises although 
the bulk phase transition is first order and not universal. The critical exponents of complete 



wetting were determined in ref . . In particular, the thickness of the confined wetting layer 
grows to a macroscopic size as 

zo(x-log(T-T,), (1.4) 
and the expectation value of the Polyakov loop at the interface vanishes as 



$1(0) oc -a/T-T,. (1.5) 



Numerical evidence for complete wetting has been reported for the SU (3) gauge theory |TD 
and for the three state Potts model [0. Complete wetting causes characteristic finite volume 
effects in the spectrum of the transfer matrix |Jl3|]. They can be used to identify complete 



wetting in a numerical simulation and to determine the value of the reduced interface tension 
(Jed- When wetting is incomplete, i.e. when add(Tc) < locd-, the transfer matrix spectrum is 
different but again allows to determine the values of the interface tensions [|1^]. Wetting is 



a phenomenon which arises only in the pure glue system. In the real world with quarks the 
Z(3) center symmetry is explicitly broken. Therefore, two of the three deconfined phases 
are only metastable. In the early universe the metastable phases have converted into the 
stable one already 10"^"^ seconds after the big bang |T^. Consequently, deconfined-deconfined 



interfaces did not survive until the QCD phase transition and wetting could not arise. To 
avoid complications due to wetting in numerical simulations of the pure glue system one 



may choose C-periodic boundary conditions ||T5|. Just like quarks they break the Z(3) 
center symmetry explicitly and eliminate two of the three deconfined phases [0. Therefore, 
in a C-periodic volume the confined phase coexists with only one deconfined phase such that 
deconfined-deconfined domain walls are absent. 



In order to apply the finite size formulae derived in ref.[13| one has to be at the finite 



volume critical temperature, i.e. at the temperature where the free energies of the confined 



and the deconfined phases are equal. In practice this temperature is not known with high 
enough precision. Therefore, in this paper we generahze the finite size formulae to temper- 
atures slightly off Tc using a dilute interface approximation. We distinguish the cases of 
complete and incomplete wetting and we also consider systems with C-periodic boundary 
conditions. For a wide class of models similar formulae have been derived rigorously by 
Borgs and Imbrie |]18| for rigid interfaces. The case of the Potts model has been worked out 



in great detail by Borgs [19|. Wang and DeTar [20| applied similar formulae to results of 



numerical simulations of the three state Potts model. 

In three dimensions interfaces living on a lattice have a roughening transition [^. For 
SU{3) pure lattice gauge theory with a fixed number of lattice points in the euclidean 
time direction we expect that the deconfined-deconfined interfaces become rigid deep in 
the deconfined phase, i.e. when the Wilson coupling (3 becomes large. A rigid interface is 
more or less fiat with only small steplike excitations following the lattice structure. Rigid 
interfaces are therefore lattice artifacts. (Of course, in condensed matter systems the lattice 
often has physical reality. Then rigid interfaces can exist in nature.) Below the roughening 
transition, which we believe to be at /3-values well above the deconfinement phase transition, 
the interfaces fiuctuate more freely. In particular, there exist soft modes — the so-called 



capillary waves — which dominate the interface fluctuations We expect that confined 



deconfined interfaces are always rough. Of course, at any given physical temperature also 
the deconfined-deconfined interfaces should become rough and hence forget about the lattice 
structure in the continuum limit. The continuum limit is also taken at large (3, but now the 
temperature remains fixed in physical units of e.g. the inverse bulk correlation length. The 
finite size formulae we are using to extract the confined-deconfined reduced interface tension 
are valid both for rough and for rigid interfaces. Still, a comparison with the results of Borgs 
and Imbrie allows to distinguish between the two possibilities. 

The paper is organized as follows. In section 2 the Z(3)-symmetry and the behavior of the 
Polyakov loop are discussed in periodic and in C-periodic volumes. In section 3 the finite size 
effects of the spectrum of the transfer matrix are derived for complete and for incomplete 
wetting as well as in the case of C-periodic boundary conditions. Section 4 contains the 
results of our numerical simulations on the subject of wetting and the determination of the 
reduced confined-deconfined interface tension acd- Finally, in section 5 we summarize our 
results and we draw some conclusions. In the appendix we give some details of the formulae 
used in section 3 to derive the spectrum in the vicinity of the phase transition. 

2 Polyakov loops in periodic and in C-periodic volumes 

Let us consider an SU (3) pure gauge theory on an x Ly x x Lt lattice corresponding to 
a temperature T = 1/Lt. The link variables U^{x,t) G 5'f/(3) are periodic in the euclidean 
time direction 

U^{x,t + Lt) = U^{x,t). (2.1) 
The order parameter for the deconfinement phase transition is the Polyakov loop 

Lt 

$(x) = Trnf^4(x,t). (2.2) 
t=x 

$(x) = $i(x) + i<^2{^) is a complex scalar field in three dimensions. The Polyakov loop is 
invariant under periodic gauge transformations g{x,t) G SU{3). However, under transfor- 
mations 

g{x,t + Lt) = g{x,t)z, (2.3) 



which are periodic only up to a center element z G Z(3) = {exp(27rm/3), = 1,2,3} it 
transforms as 

^{x)' = ^{x)z, (2.4) 

while the (Wilson) action remains invariant. Hence, the Polyakov loop indicates if the Z(3) 
center symmetry is spontaneously broken. This is the case at high temperatures. Then 
color is deconfined and an external static quark has a finite free energy F such that ($) oc 
exp(— F/T) 7^ 0. At low temperatures, because of confinement, the free energy diverges 
such that ($) = and the Z(3)-symmetry is unbroken. The two phases are separated by 
a first order phase transition At the transition temperature the confined phase coexists 
with three degenerate deconfined phases which are distinguished by the expectation value 
of the Polyakov loop. There is one phase with a real value ($) = = ($0;0). The 
other two deconfined phases are Z(3) copies of it with ($) = ^^"^^ = {—^^q, ^$o) and with 

($) = $(3) = 

In the numerical simulations we use lattices with a cylindrical geometry, i.e. with two 
short X- and ^/-directions and one much longer z-direction. Then close to the critical tem- 
perature typical configurations consist of several coexisting bulk phases, aligned along the z- 
direction, and separated by interfaces which are spanned in the transverse x- and y-directions. 
In all calculations we use C-periodic boundary conditions ITH] in the z-direction, i.e. 



U^ix + L,e,,t)= %(f,t) =[/;(x,t). (2.5) 

A C-periodic field is replaced by its charge conjugate when it is shifted over the boundary. 
Eg. (p^) implies C-periodic boundary conditions also for the Polyakov loop, i.e. 

<l>(x + L,e,) = $*(x). (2.6) 

This allows the existence of any number of deconfined-deconfined interfaces, while with 
periodic boundary conditions the number of interfaces is necessarily even. One can also use 
C-periodic boundary conditions in the x- and y-directions, i.e. 

Uf,{x + L^e^,t) = Uf,{x + Lyey,t) = U*{x,t). (2.7) 

The C-periodic action is invariant under C-periodic gauge transformations 

g{x + Liei,t) = g*{x,t). (2.8) 

Considering again a transformation which is up to a center element periodic in the euclidean 
time direction one finds 

g*{x,t) = g*{x,t + Lt)z = g{x + Liei,t + Lt)z = g{x + Liei,t)z^ = g*{x,t)z^. (2.9) 

Consistency requires = 1 and hence z = 1 (because z G Z(3)). The transformations 
of eg. ( ^.31 ) are inconsistent with C-periodic boundary conditions unless z = 1. Hence, in 
a C-periodic volume the Z(3)-symmetry is explicitly broken by the boundary conditions. 
With C-periodic boundary conditions in the x- and y-directions the Polyakov loop obeys 

$(x + L,e;) = Hx + LyCy) = $*(x). (2.10) 

In particular, the two deconfined phases and ^^^^ with non-real expectation values of the 
Polyakov loop are excluded by the boundary conditions. Hence, in a C-periodic volume the 
confined phase coexists only with the deconfined phase ^^^^ [|16|. This excludes complications 



due to wetting because deconfined-deconfined domain walls cannot even exist. On the other 
hand, when periodic boundary conditions are used in the short x- and y-directions all three 
deconfined phases coexist with the confined phase. In this situation one can study questions 
of wetting. 

In the simulations we measure correlation functions ((9j(0)(9j(z)) in the long 2;-direction. 
The operators we use are the real and imaginary parts as well as the absolute value of the 
Polyakov loop 

C»i(z)=^$i(a;,y,z), 
x,y 

0,{z)^Y.\^{x,y,z)\\ (2.11) 
x,y 

The operators Oi and are even under charge conjugation and their correlation functions 
are therefore periodic. The operator (^2, on the other hand, is C-odd and its correlation func- 
tion is hence anti-periodic. A C-odd correlation function is easier to measure in a numerical 
simulation because it gets contributions from single deconfined-deconfined interfaces while a 
nontrivial contribution to a C-even correlation function requires at least two interfaces and 
is therefore more suppressed. 



3 The interface tension and the spectrum of the trans- 
fer matrix 

The transfer matrix T{z) describes the evolution of the gluonic system over a distance z in 
the spatial z-direction. In field representation it can be defined as the transition amplitude 
between an initial Polyakov loop distribution $*^*)(a;, y) and a final distribution $*^-'^(a;, |/), i.e. 
as a path integral 

%^{z) - j VUeM-S[U]), (3.1) 

where S\U] is the gluon action. The integration is restricted to field configurations in the 
interval [0, z] for which $(a;, |/, 0) = <l'*^'-'(a;, y) and $(a;, |/, z) = ^''^\x, y). The transfer matrix 
can be written as 

r{z) ^eK^{-Hz). (3.2) 

The spectrum of H shows characteristic finite size effects which allow to distinguish between 
complete and incomplete wetting and which can be used to determine the value of the 
reduced confined-deconfined interface tension. 

For simplicity let us start with a system with C-periodic boundary conditions in the 
short X- and y-directions. Then only the deconfined phase coexists with the confined 
phase. In the subspace of these two phases the transfer matrix takes the form 

- ( til ) • 



where tdd{z) is the transition amplitude from the deconfined phase back to itself, tcd{z) 
describes transitions from the confined phase to the deconfined phase (or vice versa) and 



tcc{z) is the amplitude for transitions from the confined phase back to the confined phase. 
The eigenvalues of the transfer matrix are 

to{z) = ^[tdd{z) + t,,{z) + ^{tdd{z)-t^{z)Y + U,d{zYl 

ti{z) = ^[tdd{z) + t,,{z) - ^{Ud{z)-t,,{z)y + 4t,d{zy]. (3.4) 



In ref. |T3[ the amplitudes t^diz), tcdiz), tcc{z) were computed in the dilute interface 
approximation assuming that the free energies of all phases are equal, i.e. that one is exactly 
at the finite volume critical temperature. In practice, this temperature is not known with 
high enough precision. Therefore, to interpret the numerical data correctly, it is necessary 
to understand the transfer matrix spectrum also slightly off T^. This requires to take the 
free energy difference between the confined and the deconfined phases into account. In the 
dilute interface calculation we perform a path integral over configurations which consist of 
several bulk phases, aligned along the 2;-direction, and separated by interfaces spanned in 
the x-y-pla.ne. A block of confined phase of thickness zq costs the free energy F = fcAzo, 
where fc is the free energy density of the confined phase and A = L^Ly is the area in the 
transverse directions. Hence, such a block of confined phase is weighted with a Boltzmann 
factor 

exp(-F/T) = exp{-f,Azo/T). (3.5) 
Similarly, a block of deconfined phase gets a factor 

exp(-F/r) = expi-faAzo/T). (3.6) 
We can take out an overall factor exp (^—^{fc + fd) Azq/T^ without changing the energy 



differences. Following ref. |I9| we introduce the variable 



x = ^{f,-fd)A/T. (3.7) 

Then we write a factor exp{—xzo) for the confined phase and a factor exp{xzQ) for the 
deconfined phase. The bulk phases are separated by interfaces which cost a free energy 
F = acdAT, where acd is the reduced confined-deconfined interface tension. An interface 
gets the Boltzmann factor 

exp(-F/r) = 6exp{-a,dA), (3.8) 

where the pre-exponential factor 6 arises due to capillary wave fluctuations of the interface 
(which of course also influence the value of (Ted)- In general, one would expect that 6 also 
depends on the interface area A. However, in three dimensions the leading contribution to 
S is A-independent |T^. We treat 6 as an unknown constant. In ref.|T^ the pre-exponential 



factor was computed for rigid interfaces as S = 1. We expect that our interfaces are rough 
and that these results do not apply here. 

Let us write down the leading contributions to the amplitude tdd{z); 

tdd{z) = exp{xz) + [ dzo [ dzQ 

Jo Jzo 

exp{xzo)6 exp{-acdA) exp(-x(4 - zo))6 exp{-acdA) exp{x{z - Zq)) + ... . 

(3.9) 



The first term describes a block of deconfined phase of thickness z filling the whole interval. 
The second term is the two-interface contribution. One integrates over the positions Zq and 



z'q of the two interfaces. The first factor in the integrand stands for a deconfined bulk phase 
of thickness zq. The next factor describes the interface located at Zq followed by a block of 
confined bulk phase of thickness z'q — Zq. Next comes the factor for the interface located at z'q 
and finally a deconfined bulk phase of thickness z — z'q. Graphically one may write eq. (|3.9| ) 
as 

tddiz) d + d \ c \ d + ... , (3.10) 



and in complete analogy one constructs the expression 



+ ... . 



In fact one gets tec (-2) from t(id{z) by replacing x by — x. Similarly one writes 

/ (izoexpf— a;2;o)5exp(— o"cdv4) ex^ixiz — zq)) + ... = c 
Jo 



tcd[z) 



d + 



(3.11) 



(3.12) 



which is the one-interface contribution to the transition amplitude between the deconfined 
and the confined phase. It is straightforward, although somewhat tedious, to work out the 
multi-interface contributions and we refer the reader to the appendix for details. Using the 
variable 5' = 5exp(— o"cdA) the different terms are 



X 



X2 + 

X 

x2 + 6'' 



f_ 

x2 + 5'^ 



sinh(zY x"^ + 5' 



tddiz) = cosh{z\/ x^ + 6'') + 
tcc{z) = cosh(z-\/ + 5'^) — 

tcd{z) = 

leading to the eigenvalues 

to{z) = exp(z\/a;2 + 6'^), 
ti{z) = exp{—z\/x^~+~S^). 

According to eq.(|3.2|) we define Ei by 

ti{z) = exp{-Eiz), 

and we find for the energy difference 



smh.{zyx'^ + 6' ), 



smh.{zy x"^ + S' ), 



(3.13) 



Ei-Eq = 2Jx^ + 6^exp{-2acdA). 



(3.14) 



(3.15) 



(3.16) 



At the finite volume critical temperature (i.e. for x = 0) the result agrees with ref.[|l^]. It 
predicts that the energy splitting between the lowest states has a minimum at the critical 
temperature. At the minimum the energy splitting 



Ei-Eq = 25exp(-(TedA) (forx = 0) 



(3.17) 



vanishes exponentially with the interface area times the reduced confined-deconfined inter- 
face tension acd- This equation will be used to extract the value of acd in the numerical 
simulations. 



Now let us turn to the more complicated case of periodic boundary conditions in which 
the confined phase coexists with all three deconfined phases. As it was shown in refs. 19 1 
in the subspace of the three deconfined phases ^^^^ and the confined phase the 

transfer matrix takes the form 



t{z) 



( tdd{z) tdd'{z) 

tdd'{z) tdd{z) 

tdd'{z) tdd'{z) 

V tcd{z) tcd{z) 



tdd'{z) tcd{z) \ 

tdd'{z) t^d{z) 

tdd{z) t^d{z) 

tcd{z) tcc{z) J 



(3.18) 



where tdd{z), tdd'{z), tcd{z) and tcc{z) are transition amplitudes between the various phases. 
The transfer matrix t{z) is symmetric under the permutation group = Z(3) x C which is 
only broken by the boundary conditions in the z-direction. This group has three irreducible 
representations, namely a symmetric (trivial), an antisymmetric, and a mixed symmetric 
one, corresponding to the Young tableaux 



^ (3.19) 

It turns out that the spectrum of the above transfer matrix splits up into a symmetric ground 
state with the eigenvalue to{z), a corresponding excitation t^i^z), and a twofold degenerate 
eigenvalue ti,2(-2) of mixed symmetry. The eigenvectors are of the form 

Vs = {vd, Vd, Vd, Vc) (3.20) 
for the symmetric representation, and 

3 

Vm = {vi,V2,V3,0) with "^Vi = (3.21) 

i=l 

for the mixed representation. For the latter we chose one basis vector |1) with positive C- 
parity and the other one |2) with negative C-parity. The corresponding eigenvalues of the 
transfer matrix are 

to{z) = ^[tdd{z) + 2tdd'{z) + tcc{z) + ^{tdd{z) + 2tdd'{z) - t,c{z)f + l2t,d{zYl 
^1,2(2;) = tdd{z) — tdd'{z), 

h{z) = ^[tdd{z) + 2tdd'iz) + tcc{z) - ^{tdd{z) + 2tdd'{z) - tUz)y + I2t,d{zy]. 

(3.22) 

Due to the S's-symmetry the eigenvalues ^1(2;) and t2{z) are degenerate. We must distinguish 
two cases: complete and incomplete wetting. 

Let us begin with the complete wetting case. Then a deconfined-deconfined domain 
wall consists of two widely separated confined-deconfined interfaces. Therefore, in the dilute 
interface approximation one must include only these interfaces, while direct deconfined- 
deconfined interfaces do not exist. In graphical notation the various amplitudes read 



tdd{z) 
tdd'{z) 

tcd{z) 
tcc{z) 



+ 



d 



d + ... 



d 



d^ + 
+ ... , 

+ 3~ 



+ ... . 



(3.23) 



Due to complete wetting the leading contribution to the transition amplitude tdd'{z) between 
two different deconfined phases comes from two confined- deconfined interfaces. The factor 3 
in the amplitude tcc(^) arises because there are three possible realizations of the deconfined 
phase between the two blocks of confined phase. Following the appendix we now put 5' = 
■\/3(5exp(— (TcdA). Working out the multi-interface contributions results in 



V^tcd 




5' 








tec 




= cos\i{z\J X 


tdd{z) + 2tdd' 




= cosh(z\/ X 


tdd{z) — tdd' 




= exp{zx), 



sinh(2;\/ + 6' 



2\ 

5 



2 + 5'2 



; sinh(2;Y x^ + 5' 



x2 + 5'^ 



sinh(z\/ x^ + 5' 



(3.24) 

where tcc{z) and tdd{z) + 2tdd'{z) are related by exchanging x with — x. The eigenvalues are 



to{z) = exp(z\/x2 + 5'^), 
^1,2(2;) = exp(zx), 



t3{z) = exp(-2\/x2 + 5'2). (3.25) 



resulting in the energies 



^1,2 - ^0 = V^' + 352exp(-2aedA) - x, 



E3-Eo = 2a/x2 + 352exp(-2acdA). (3.26) 



At the critical temperature (for x = 0) the results are in agreement with ref.|13]. To leading 



order in x they are identical with the finite size formulae derived by Borgs [T^ who neglected 
higher-order multi-instanton contributions. We note that the x-dependence disappears in 
the combination 

{E3 - ^1,2) (^1,2 - ^0) = 35^exp{-2a,dA). (3.27) 



Eq. (|3.26|) predicts that the energy E^—Eq has a minimum at the critical temperature (x = 0), 



at which 

£^1,2 - Eq = exp(-crcd^), 

- = 2V35exp(-(TcdA) (forx = 0), (3.28) 

such that one can again read off the value of acd- The equidistant energy splitting between 
the lowest levels at the critical temperature, i.e. 

2(Ei,2 - ^0) = ^3 - ^0 (forx = 0), (3.29) 

is characteristic for complete wetting. It disappears in the incomplete wetting case to which 
we now turn. 

In case of incomplete wetting also direct deconfined-deconfined domain walls exist, 
which do not consist of two confined-deconfined interfaces. They have another reduced 
interface tension add and they get a Boltzmann factor 



exp(-F/r) = 7exp(-arfrfA). 



(3.30) 



The various amplitudes then get additional contributions 



tdd'(z) 

tcd{z) 



= d + d 


c 


d + 2 d 


d' 


d + 




= d 


d' + d 


c 


d' + d 


d" 


d' + 






— c 


d + 2 c 


d' 


d + ... , 




= c + 3 c 


d c + ... . 



(3.31) 

In particular, the leading contribution to the amplitude tdd'{z) now comes from a direct 
deconfined-deconfined interface. Summing up all multi-interface contributions and using the 
variables 7' = 27exp(— 0"^^^) and 5' — \/2>5 exp(— (Jcd^) along with x' — x -\- 7'/2 yields 



8' 



smh.[z\l x'"^ + 5'^) exp(2;7'/2), 



tcc.iz) 



cosh(^Yx'^ + 5'"^) — 
tddiz)+2tdd'iz) = \cosh{z\/x'^ + S'^) + 



X' 



x'^ + 6' 



sinh(^V2;'^ + ^ ) exp(^772), 



X 



+ S'' 



smh{z\/ x'^ + 5'^) I exp(z772), 



tdd{z) — tdd'{z) 



exp(^a;) exp(— ^7'/2). 



(3.32) 



We arrive at the eigenvalues 



^0(2;) = cxp{z\/ + 5'^) cxp(2;7'/2), 
^1,2(2;) = exp(2;a;) cxp(-2;772), 

tsiz) = exp{-z\/x'^ + 6'^) exp{z'y'/2) 



(3.33) 



and hence 



-^1,2 - -Bo = \/ (x + 7exp(-(7dd^))2 + 352 exp(-2(7cd^) - x + 27exp(-(7ddA), 
E3-Eo^ 2^{x + 7exp(-(7dd^))2 + 35^ exp(-2(7edA). (3.34) 

The combination 

(^3 _ Ei,2)(^i,2 - ^0) = 35' exp(-2aedA) - 37' exp(-2arfdA) + 60:7 exp(-add^) (3.35) 

is now x-dependent. Furthermore, the minimum of E3 — Eq no longer occurs at the critical 
temperature. Instead it corresponds to 

X = -^exp{-addA) < 0, (3.36) 

which lies in the confined phase, because x < implies that the free energy of the confined 
phase is smaller than the one of the deconfined phase. At the minimum one finds 

Ei,2 -Eo^ VSS exp{-acdA) + 37exp(-(7<idA), 

E3-Eo^ 2V35 exp(-(7edA) (for X = -7 &q){-addA)). (3.37) 

As opposed to the complete wetting case the point of equidistant energy splittings between 
the lowest levels is now at a different position 



X = 2'y exp{-(TddA) > 0, 



(3.38) 



which lies in the deconfined phase. It is clear that the transfer matrix spectrum is qualita- 
tively different from the complete wetting case. 



4 Numerical results 



We have performed numerical simulations of the SU{3) pure gauge theory on lattices with 
Lt = 2 points in the euclidean time direction close to the deconfinement phase transition. 
For the standard Wilson action the critical coupling is then given hj (3 = 5.0933(7) p3| . 
The spatial geometry is cylindrical. We have used the sizes x Ly x Lz = 4: x A x 64, 
4 X 6 X 64, 6 X 6 X 64, 6 X 8 X 96 and 8 x 8 x 128. In all cases C-periodic boundary 
conditions have been used in the long z-direction. In the transverse x- and y-directions we 
have used both periodic and C-periodic boundary conditions. Of course, thermodynamics 
requires periodic boundary conditions for the gluon field in the euclidean time direction. We 
have performed measurements at several /3-values very close to the critical coupling using 
an overrelaxed heat bath algorithm. To interpolate our results to other /3-values we have 
used standard reweighting techniques. Per lattice we have typically performed 50000-250000 
measurements. 

The low energy states in the spectrum of the transfer matrix in the 2;-direction are 
extracted from the correlation functions of the operators Oi of eq.( |2.11D . Since we always 
use C-periodic boundary conditions in the ^-direction we introduce a complete set of states 
\n) with their charge conjugate states C\n) and we obtain 

Z{O,{0)O,{z)) = Y.(^\0,exp{-Hz)0,exp{-H{Lz-z))C\n) 

n 

= ^{n\Oi\m) exp{-E^z){m\Oi\n) exp {-En{Lz - z)) Cn 

m,n 

= ^^CnlinlOilm)]^ exp{-Emz) exp {-En{Lz - z)) , 

m,n 

Z = Y.(^\exv{-HLz)C\n) = Y.Cnexv{-E^Lz), (4.1) 

n n 

where we have used C\n) = Cn\n) for states of C-parity c„ = ±1. 

For C-periodic boundary conditions in the x- and y-directions only the operators Oi and 
O3 are relevant. Both create the state with energy Ei. Putting Eq = and denoting the 
states with energies Eq and Ei by |0) and |1) the transfer matrix formalism predicts 

Z{Oi{0)Oi{z)) ~ mOi\l)\\expi-E,z) + exp{-E,{L,-z))) 
+ |(0|Oi|0)p + |(l|Oi|l)|2exp(-^iL,), 

Z{Os{0)O,{z)) ~ mOs\l)\\expi-E,z) + expi-E,iL,-z))) 
+ K0|O3|0)P + Kl|O3|l)|'exp(-EiL,), 
Z ~ l+exp(-EiL,), (4.2) 

such that 

{O,{0)Oi{z)) ~ Acosh{E,{z-hz)) + B, 

{O,i0)Osiz)) ~ Ccosh{E,iz-hz)) + D. (4.3) 

We use these equations to determine the energy Ei by fitting the measured correlation 
functions. 

When periodic boundary conditions are used in the x- and y-directions the energy split- 
tings Ei^2 — Eq and E^ — Eq between the excited states |1), |2) and |3) and the ground state 



|0) become small close to the phase transition. Applying the reduction rule 



n 



(4.4) 



for 5*3 (see [22|) and again setting Eq = we obtain 



+ mo, 

X (exp(- 

Z{02{0)02{Z)) ~ |(0|O2 

- 

X (exp(- 

Z{Os{0)O3{z)) ~ |(0|O3 

+ |(0|O3 



1) \\exp{-E,z) + exp(-Ei(L, - z))) 
3)pexp(-EiL,) 

■(^3 - E,)z) + exp(-(E3 - Ei)(L, - z))), 

2) \\exp{-E2z) - exp(-E2(L, - z))) 

3) pexp(-E2L,) 

-(^3 - E2)z) - exp(-(E3 - E2)iL, - z))), 
3)p(exp(-E3z) + exp(-E3(L, - z))) 
0)p + |(3|O3|3)pexp(-E3L,), 



1 + exp(-E3L^ 



(4.5) 



Here the additional symmetry relation | p = | (2|(9j|2) p was used. Note that the corre- 

lation function of the C-odd operator O2 is anti-periodic because of C-periodic boundary con- 
ditions in the z-direction. The S'3-symmetry ensures that E, = E2, |(0|(9i|l)p = |(0|(92|2)P 
and |(l|Ci|3)p = | (2|C»2|3) p. For very large L, such that EiL, > 1 eq.(|3D simplifies to 



m0)O,{z)) 
{O2i0)O2iz)) 
mO)Os{z)) 



Acosh{Ei{z-^L,)), 
Bsmh{E2iz-h,)), 
Ccosh{E,{z-h,)) + D. 



(4.6) 



In most of our fits we have used this simpler form. However, in some cases we have checked 
our results using the more general expressions ( |4.5| ) and we found agreement within the error 
bars. 

Fig.l shows a configuration on the 8 x 8 x 128 x 2 lattice which consists of several 
bulk phases separated by interfaces. Two typical correlation functions are depicted in fig. 2. 
We should mention that autocorrelation times of the correlation functions are rather long, 
especially at large distances. Therefore, it was always a good check if the energies E, and 
E2 agreed within errors (as they should because of 5'3-symmetry) . 

Fig. 3a shows the energy i?i on a 4 x 6 x 64 x 2 lattice with C-periodic boundary conditions 
in the x- and y-directions. As expected, there is a minimum at the critical coupling. From 
the finite size effect of the value at the minimum we have determined the reduced confined- 
deconfined interface tension using eq.( p.l7|) . The result of the fit, which is depicted in fig. 4, 
is given by 

acd = 0.035(1) = 0.140(4)T2 (C-periodic). (4.7) 



The fitted value of the pre-exponential factor is 6 
prediction 6 = 1 for rigid interfaces 
deconfined interface is indeed rough. 



0.196(5) which is much smaller than the 
We take this as an indication that the confined- 



Figure 1: The z-dependence of Oi{z) = I]^,^ $i(a;, 2) for a configuration on an 
8x8x128x2 lattice with C-periodic boundary conditions in the spatial directions at = 
5.0933. 



Figure 2: The correlation functions (Ci(O)Oi(z)) and (02(0)02(2)) ona4x6x64x2 
lattice with C-periodic boundary condition only in the 2;-direction. The solid lines represent 
fits of the data to the corresponding equation (|4.6|) . 



Figure 3: The spectra of the lowest energies in the region of the phase transition tem- 
perature. In a) the energy Ei from simulations ona4x6x64x2 lattice with C-periodic 
boundary conditions in the x- and y-direction is shown. In b) the degenerate energy E12 and 
the energy derived from simulations ona6x6x64x2 lattice with periodic boundary 
conditions in the x- and y-direction are shown. The energies are plotted as boxes at the 
simulation points, whereas the interpolated values are plotted as diamonds. 



Figure 4: The dependence of the energies on the cross section A = L^Ly of the cyhndrical 
lattices. The values of log(i?i) are shown from simulations on lattices with C-periodic 
boundary conditions in the x- and y-direction (diamonds) and the values of log(ii^i,2) (= 
log(^i?3) ) are shown from simulations on lattices with periodic boundary conditions in the 
X- and y-direction (squares). The solid lines represent the fits of the data to eg. ( P.17| ) and 
eq.( |3.28| ) respectively. 



Fig.3b shows the energies £'1,2 and ii^a on a 6 x 6 x 64 x 2 lattice with periodic boundary 
conditions in the x- and y-directions. As expected, the energy has a minimum which is, 
however, shghtly shifted away from the critical (3 into the deconfined phase. At the same 
time, the point where 2Ei 2 = is shifted slightly into the confined phase. In case of 
incomplete wetting one would expect such shifts, but they would go exactly in the opposite 
direction. Therefore, our data favour complete wetting. Of course, then eq.( |3.26| ) predicts 



that the minimum of £"3 and the point where 2i?i 2 = -E3 should correspond to the same 
/3-value, which is almost — but not quite — the case for our data. We believe that this 
effect is some sub-leading finite size effect which is not correctly reproduced by our dilute 
interface calculation. The same effect is visible in the data of Wang and DeTar in the 



three-dimensional three state Potts model, which also shows complete wetting There 
the shift is present on smaller lattices, and it disappears on larger volumes, supporting the 
assumption of a sub-leading effect. Altogether, we consider fig. 3b as a strong indication 
that the deconfined-deconfined interfaces are completely wet by the confined phase at the 
deconfinement phase transition. To extract the reduced confined-deconfined interface tension 
from our periodic lattices we have taken the value of Ei^2 at the point where 2Ei^2 = E3 
and we have fitted its finite size behavior to eq. ( ^.281 ). The fit, which is also shown in fig. 4, 
results in 

acd = 0.034(2) = 0.136(8)T2 (periodic). (4.8) 

It is reassuring (and further supports the assumption of complete wetting) that this result 
is in agreement with the one obtained from the C-periodic simulation. Also the fitted 
value of the pre-exponential factor, S = 0.202(6), agrees with the one obtained from the 
C-periodic data. This value again indicates that the confined-deconfined interface is rough. 
Furthermore, it means that the capillary waves have the same behavior in periodic and in C- 
periodic volumes. This was expected because they are small fluctuations which are insensitive 
to global topological characteristics of the boundary conditions. Taking the results of both 
boundary conditions together we state as our final answer for the reduced confined-deconfined 
interface tension 

= 0.139(4)Tf . (4.9) 

We have also re-analyzed the data of Wang and DeTar from our perspective. They 
have measured the energies Ei^2 and E^ in the three-dimensional three state Potts model 
using periodic boundary conditions. As Wang and DeTar have pointed out, their results are 
in good agreement with complete wetting. On a 20^ x 120 lattice they find £^12 = 0.048(1) 
at /3 = 0.3668 and ^1,2 = 0.0397(6) at (3 = 0.3669 while at = 0.3670 on a 30^ x 120 
lattice Ei^2 = 0.0185(3). In all cases 2Ei^2 = E3 within the error bars. Assuming that the 
asymptotic formulae (|3.28| ) are valid for these lattices we obtain 

(Ted = 0.0017(3) (Potts model), (4.10) 

where the error is mostly systematical. This is due to the same slight shift of the minimum 
of £"3 away from the point where 2£i 2 = £3 which we also observed in our data. In the 
Potts model the shift is present only on the smaller lattice, while it disappears on the larger 
volume. Perhaps the result of eq.( 4.10| ) should be taken with a grain of salt because it relies 



only on two different volumes. For the pre-exponential factor we find 6 = 0.042(4) indicating 
that the ordered-disordered interface is rough. The ordered-ordered reduced interface tension 
Cdd{Tc) has been determined in the Potts model by Karsch and Patkos |jTl|. They enforced 



the interface by boundary conditions and computed the free energy difference due to the 
interface by a numerical integration. Their result 



cTdd{Tc) > 0.0026(2) (Potts model) 



(4.11) 



(translated to our notation) is consistent with eg. ( [4. 101) because complete wetting implies 



5 Conclusions 

We have related the finite size effects of the spectrum of the transfer matrix to the value 
of the reduced confined-deconfined interface tension both for C-periodic and for periodic 
boundary conditions. In the periodic case one must distinguish between the possibilities of 
complete and incomplete wetting. Our numerical data strongly suggest that the deconfined- 
deconfined domain walls of the SU (3) pure gauge theory are completely wet by the confined 
phase. Hence, in a hypothetical universe, which contains only gluons but no quarks (or 
other particles), the deconfinement phase transition would proceed via domain wall splitting 
and not via the usual bubble nucleation. However, the presence of quarks in our universe 
made wetting at the deconfinement phase transition impossible |]12|, [1^. Still, the reduced 
confined-deconfined interface tension of the pure glue system serves as an upper limit on the 
corresponding value in real QCD. The determination of the value ad = 0.139(4)T^^ required 
a careful analysis of finite size effects and an understanding of the wetting dynamics. 

We like to mention that the transfer matrix method is a very powerful tool to extract 
interface properties from numerical simulations, especially because finite size effects are 
well understood. Certainly, in this method information on correlation functions at large 
distances is needed. With only local algorithms available (as it the case for gauge theories) 
this is not an easy task, and it requires very high statistics. An attractive alternative is the 



multicanonical algorithm of Berg and Neuhaus p4|. We have also applied this algorithm to 



the pure SU{3) gauge theory |^ and we have measured the interfacial free energy. It has 
turned out that the use of cylindrical lattices is again advantageous in order to reduce the 
interactions between interfaces. This work is presently in progress. 
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Appendix 

In this appendix we would like to give some of the details for the formulae occurring in 
section 3. We shall start with the complete wetting case. Let us also assume that we 
have q deconfined phases, such that g = 3 for periodic and g = 1 for C-periodic boundary 
conditions. It is convenient to define the following integrals: 

= diJ di2--- di^e^zyY^^-im n>l (A.l) 



where In{y, 1) = and y = 2x = {fc — fd)A/T. The aim is to evaluate /„(y) = I{y, 0), which 
exactly corresponds to n confined-deconfined interfaces. For reasons of notational simplicity 



we rescale zx x, zy ^ y, and z^yq6 exp{—acdA) = z6' — > 6'. The different interface 
contributions ^/qtcd, tec and tdd + {q — ^)tdd' can now be expressed in terms of In{y)'- 



tec 

tdd + {q- 

The case n = 1 is trivial and yields 



exp(a;)^/2.+i(y)5'^"+S 

n=0 

oo 

exp(-x)(l + ^/2n(z/)0, 

n=l 
oo 

exp(x)(l + ^72n(-y)0- 



n=l 



1 



h{y,0 = -(exp(-yO -exp(-y)), 
= -(1 -exp(-y)). 

y 

For n > 2 we find after differentiating /^(y, ^) with respect to ^ 



In particular for n = 2: 



The integration is easy and we obtain 



-exp(-y04-i(-y,0- 



-(l-exp[y(l-0]). 

y 



l2(y,0 = ^ + ^(exp[y(l-0]-l), 
Hy) = -i + :^j(exp(y)-l) = -exp(y)-^^^^^^ 



y y 

Repeating the steps above we arrive at 



dy 



hn{y) = exp(|/) , 

n dy 

hn+\\y) = -exp(-y)— , n>l. 

n ay 

This we checked explicitly up to n = 3. Prom the previous two equations follows: 



hn+i{y) = 

l2n+2{y) = 



(-1)" 

(-1) 



dy 



n!(n + l)! 



d ^2 

-7^ + 



9y dy^ 



hiy), n>l. 



Now we can do the -^/qtcd term. Using the identities: 



4" " 

exp(a;)/2n+i(2a;) = -r^Yl 



i=0 



n 
i 



.Ni <9^* / / 1 5" /sinhx' 
-1 exp a; /i 2a; = -tttt^ 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



we find 

V^tc, = E f ^) S"^'^' = sinh V^^. (A.IO) 

Now we come to tec- Using eqs. (|A.2|) and ( |A.7] ) we can write: 

tec = exp(-a;) - exp(x)— ^ — /2„_i(y)5'2". (A.ll) 
Differentiating the previous equation with respect to 6' and integrating we obtain 

tec = exp(-x) - exp(x) d6'—{exp{-x)^/qtcd) 

= + , ^ - sinhVx^ + 6'^. (A. 12) 

Likewise we get 



tdd + {q- 1)W = cosh Vx^ + 5'2 + sinh + 6'^. (A. 13) 

Finally we have trivially 

tdd - tdd' = exp(x). (A. 14) 

The case of incomplete wetting is somewhat harder. First we rescale the expression 
z{q — l)7exp(— cTrff^A) = 27' — 7'. We begin again with ^tcd- To lowest non-trivial order 
in 5', but to all orders in 7', we obtain easily 



, • n « 
^1=0 



(A.15) 
The single sum is easy 

For the double sum we first let i — > n—i and then we interchange the order of the summations. 
This gives 

ig/ yyv^ (-i/r^ _ ly^/ yy^ _ exp(y) ^^^^^ 

so that with x' = x + 7'/2 

v^tcd = 5'^^exp(y/2). (A.18) 
Combining eqs. ( |A.1(]| ) and ( [A.17| ) we conjecture that to all orders in 5' the following holds: 

5' 



^^'^ = TPtF "'""^ v^^" + '^"exp(772). (A.19) 

This we checked explicitly up to order 0{6'^). It is now easy to guess the expressions for t^c 
and tdd + {q- ^)tdd'- They are: 

tec = (coshVx'^ + 6'^ - ^ ^ - sinhVx^2^^^2[exp(772), 



tdd + {q- 1)W = ( cosh + 5'2 + sinh + 6'^ ] exp(y/2). (A.20) 



This we also checked explicitly up to order 0{5'*). Notice for 7' = we reproduce the results 
for the complete wetting case. Finally we find 

tdd - tdd' = exp(a;) exp[-77(g - 1)]. (A.21) 

The formulae given in the main text are obtained by replacing 

X — >• zx, x' — zx\ 5' z5' = ^^5cxp(— (Tcd^), 7' — z'^' = z{q—l)'j exp{—acidA). (A. 22) 

As an addendum we would like to mention the following identity: 

det(i) = exp[-(/, + qfd)A/T], (A.23) 

which is equivalent to 

tcc{tdd +{q- - ^/qtcd^/qtcd = exp(7'). (A.24) 
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